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Abstract. Low-luminosity active galactic nuclei (LLAGNs) represent the bulk of the AGN 
population in the present-day universe and they trace low-level accreting supermassive black 
holes. The observational properties of LLAGNs suggest that their central engines are intrinsically 
different from those of more luminous AGNs. It has been suggested that accretion in LLAGNs 
occurs via an advection-dominated accretion flow (ADAF) associated with strong jets. In order 
to probe the accretion physics in LLAGNs as a class, we model the multiwavelength spectral 
energy distributions (SEDs) of 24 LINERs (taken from a recent compilation by Eracleous et al.) 
with a coupled accretion-jet model. The accretion flow is modeled as an inner ADAF outside of 
which there is a truncated standard thin disk. These SEDs include radio, near-IR to near-UV 
HST data, and Chandra X-ray data. We find that the radio emission is severely underpredicted 
by ADAF models but can be explained by the relativistic jet. The origin of the X-ray radiation 
in most sources can be explained by three distinct scenarios: the X-rays can be dominated by 
emission from the ADAF, the jet, or from both components contributing at similar levels. From 
the model fits, we estimate important parameters of the central engine of LINERs, such as the 
mass accretion rate - relevant for studies of the feeding of AGNs - and the mass-loss rate in the 
jet and the jet power - relevant for studies of the kinetic feedback from jets. 
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1. Introduction 

Th anks in large part to the Palomar survey (e.g., iHo. Filippenko fc Sargent 1995 



19971) . we know that most active galactic nuclei (AGN) in the present-day universe 
have low luminosities, being thus called low-luminosity AGNs (LLAGNs). The bulk of 
the LLAGN population (w 2/3) are low-ionization nuclear emission-line regions (LIN- 
ERs), which are e xtremely sub-Eddington systems with an av erage Eddington ratio of 
iboi/iEdd ~ 10" 5 (|Holl2009l: lEracleous. Hwang fc Flohidl201fi . 



The observational properties of LINERs - and LLAGNs in general - are quite different 
from those of more luminous A GNs. Regarding the SEDs , LLAGNs seem not to have 
the big blue bump feature (e.g.. iNemmen et all [20(51 IHoI l2008t lYuan, Yu fc Hoi 120091 : 



but see lMaojl2007h which is one of the signatures of the presence of an optically thick, 



geometrically thin accretion disk. Regarding the emission -lines, LLAGNs typically have 



weak and narrow Fe Kcv emission ( Terashima et al.l 2002 ) and a handful of LINERs dis 



play broad double-peaked Ha lines (e.g., Storchi-Bergmann et al.ll2003l ): these properties 



of the emission-line spectrum are consistent with the absence of a thin accretion disk, 



119 



120 



R. S. Nemmen, T. Storchi-Bergmann, M. Eracleous & F. Yuan 



or a thin accretion disk whose inner radius is truncated at > 100GM/c 2 . Last but not 
least, with the typical fuel supply of hot diffuse gas (via Bondi accretion) and cold dense 
gas (via stellar mass loss) available in nearby galaxies, LLAGNs would be expected to 
produce much higher luminosities tha n observe d on the assumption of standard thin 
disks with a 10% radiative efficiency (Ho 2009() . Taken together, this set of observa- 
tional properties favors the scenario in which the accretion flow in LLAGNs/LINERs is 
advection-dominated or radiatively inefficient. 

Ad vection-dominated accretion flows (ADAFs; for a recent review see lNaravan fe McClintock 



20081) are very hot, geometrically thick, optically thin flows which are typified by low ra- 
diative efficiencies (L <C O.lMc 2 ) and occur at low accretion rates (M < O.OIA/Edd)- 
Supermassive black holes are thought to spend > 95% of their lifes in the ADAF state 
( Hopkins. Naravan fc Hernquistl 12006 ). the best studied case being Sgr A* (e.g., lYuanl 
20071) . 

ADAFs are relevant to the understanding of AGN feedback since they are quite efficient 
at producing powerful outflows and jets, as sug gested by theoretical studies (including an- 



alytic al theory and numerical simulations; e.g.. lNemmen et al.ll2p07|lNarayan fc Mc Clintock 



2008), as well as different observational studies of LLAGNs (e.g., Hoi2002t Heinz. Merloni fc Schwab! 
2003). ln fact, the so-called "radio mode" of AGN feed back invoked in sem i-analytic and 
hydrodynamic simulations of galaxy formation (e.g., Bower et al.l l2006h would corre- 
spond to the A DAF accretion sta t e actively producing jets as explic itly incorporated 
in some works ( Siiacki et al. 120071 : lOkamoto. Nemmen fc Bower 2008). An alternative 
and perhaps more appropriate expression for the "radio" feedback mode would be the 
"LINER/ ADAF" mode. 

It is clear that advances in the understanding of the physical nature of LINERs and 
LLAGNs are required in order to understand the nature of black hole accretion and 
feedback in the local universe. The goal of this work is therefore to probe the physics 
of accretion and ejection in the LINER population, by modeling their nuclear multi- 
wavelength SEDs which provide constraints to physical models for the emission of the 
accretion flow and the jet. 



2. Data set 

Our data set consists of 24 SEDs of LINERs which include radio (VLA), near-IR - 
optical - UV (HST) and X-ray (Chandra) data with high spatial resolution. These SEDs 
were selected from a sample of 35 SEDs compiled by|| racleous. Hwang fc Flohid (<2010h 



using two selection criteria: (i) there should be estimates of the black hole mass for the 
corresponding galaxies and (ii) there should be good X-ray estimates of the photon index 
and X-ray luminosity. 

Based on these criteria, we can separate the SEDs in two groups: group A, comprising 
10 LINERs with the most complete sampling of the SEDs, and group B, for which there 
is lack of data in some parts of the SED. For illustration! we list the LINERs in the 
group A: NGC 1097, M81, NGC 3998, NGC 4143, NGC 4278, M84, M87, NGC 4579, 
NGC 4594 and NGC 4736. 



3. Models for the continuum emission 

In order to model the LINER SEDs and constrain the properties of their central en- 
gines, we adopt th e physical scenario wh ich is favored to explain t he observational prop- 
erties of LLAGNs (|Yuanll2007t lHoll2008t see Figure 13 of [Holl2008l for a cartoon). In this 
model, the accretion/ejection flow consists of three components: (i) the inner parts of the 
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flow are advection- dominated, geometrically thick; (ii) the outer parts of the accretion 
flow are in the form of a standard thin disk truncated at a certain transition radius; (iii) 
near the innermost parts of th e ADAF a relativistic jet is launched. More details can be 
found in Yuan. Yu &: Ho 20091 and Nemmen et al. 2010 (in preparation). 

The radiative processes operating in the ADAF are synchrotron emission, bremsstrahlung 
and inverse Compton scattering. The truncated thin disk spectrum is simply thermal. 
The jet contributes with synchrotron emission. Some of the main parameters of the 
models will be discussed in the section below. 



4. Results 

We describe the results of the spectral fits for only two LINERs of our data set for 
the sake of brevity: NGC 4374 (M84) and NGC 4594 (Sombrero). These two examples 
nevertheless are illustrative of our general results. In the SED plots below, the arrows 
in the IR correspond to upper limits due to significant contamination of the nuclear 
emission by the host galaxy. The error bars in the near-IR to UV denote the uncertainty 
due to the range of possible extinction corrections. 

When fitting the SEDs, we explored the parameter space of the accretion/jet models 
within the range of values plausibly allowed by theory. We avoided to fix the values of 
some dynamical and microphysical parameters for which there is a substantial uncertainty 
on theoretical grounds. One example is the parameter 5, which controls the amount 
of energy dissipated via turbulence that is directly deposited on electrons, for which 
there is a considerable u ncertainty with plausible values in the range ~ 0.001 — 0.5 



(e.g., ISharma et al.l 120071) . In all models we considered that only some fraction of the 



mass accretion rate available at the outer boundary of the accr etion flow ends up bein; 



accre ted, due to mass-loss via winds produced in the ADAF (jNaravan fc McClintoc 

[200I). 

The SED of M84 is plotted in Figure Q] together with two different spectral fits. The 
left panel of this figure shows a model in which the ADAF dominates the emission from 
the near-IR upwards, particularly the X-ray emission. The thin disk dominates the flux 
output around ~ 10/mi and the jet is responsible for the bulk of the core radio emission. 
For this model, the mass accretion rate supplied at the outer radius R tI of the ADAF 
is M(R tr ) = 4 x 10~ 4 A/Edd, compatible with the Bondi accretion rate estimated by 



Pellegrinil ([20051 ); furthermore, R tr — 150i?s and S = 0.01. The mass-loss rate in the jet 
was estimated as Mj et = 10 _3 M(i? t r)- 

The right panel of Figure Q] shows a model for the SED of M84 in which the jet 
dominates the radio and X-ray emission. Notice that the shape of the X-ray spectrum 
from the ADAF is not consistent with the data. For this model, M(R tI ) = 1.5 X 10~ 4 MEdd 
(again compatible with the Bondi accretion rate), R tI — 30i?s, 8 — 0.3 and Mj ct = 

0.01M(.Rt r ). 

Therefore, we demonstrated that there are two possible types of models which can 
accomodate the observed SED of M84. In the first type, the emission from the ADAF 
dominates the observed X-rays; in the second type of model, the jet emission dominates 
the X-rays. By construction, a third type of model is also possible in which the jet and 
the ADAF contribute with similar intensities to the high energy emission. These results 
apply also to the other LINERs in our sample. To illustrate this, we also show spectral 
fits to the nuclear SED of the Sombrero galaxy in Fig.[5J The left panel shows an "ADAF 
dominated X-rays" model with M{R tr ) = 9.1 x 10" 4 M E dd, Rtr = 300R S , S = 0.01 and 
Mjet — W M(Rtt)- By varying the microphysical parameters of the jet model, we are 
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Figure 1. Models for the SED of M84, including the contribution of the ADAF, truncated thin 
disk and relativistic jet. Left panel: model in which the ADAF dominates the observed X-ray 
emission. Right panel: model in which the jet dominates the X-ray output. 
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Figure 2. Same as Figure[T]for NGC 4594 (Sombrero). Only a jet model is shown in the right 
panel, for which the mass- loss rate is Mj e t = 4.5 x 10~ 7 MEdd- 



also able to obtain a "jet-dominated X-rays" model which also successfully explains the 
entire SED (right panel). 

Figure El shows the qualitative aver age SED of LINERs based on the one computed 
by i racleous. Hwang fc Flohid (2010). From our modelling of the SEDs of LINERs, 
we are able to unveil the physical nature of the LINER continuum emission in each 
waveband, as outlined in the upper part of Fig. [3] The uncertainty about the ori- 
gin of X-rays in LINERs has been debated in the context of LLAGNs and also Sgr 
A* (Merloni. Heinz fc Di Matteoll2003t iFalcke. Kordinefc" Markofj[2003lYuan. Yu fc Hoi 
2003 ) . In our case, the root of the uncertainty in the origin of X-rays lies in the uncer- 



tainties regarding the microphysics of the hot plasma in the ADAF and the jet (mainly 
the uncertainty on the value of 5 and the effect of shocks in the jet), which has a ma- 
jor impact in the fitting of the observed X-ray spectrum. We suggest that monitoring 
campaigns of the variability of radio and X-ray emission in LINERs, and the comparison 
of such observations with predictions of jet/ ADAF models would help to pin down the 
nature of the X-ray emission in LLAGNs. 

From our modelling of the LINER SEDs with the coupled jet-ADAF model, we are 
able to constrain important parameters that characterize their central engines: 

• The typical accretion rate available at the outerskirts of the accretion flow is M ~ 
10- 3±1 M Edd . 
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Figur e 3. Qualitative average SED of LINERs (solid line) based on Eraclcou s. Hwang fc Fl ohic 
l|2010D . The dashed line illustrates the potential effect of extinction corrections. The top part 
of the plot shows which components of the flow dominate the emission in each waveband. Also 
displayed are the averag e SEDs of radio-loud (dotted line) and radio-quiet (dot-dashed line) 
quasars from Elvis (1994) for the sake of comparison. 



• Typical jet mass-loss rates are in the range Mj ot ~ 10 8 — 10 4 AfEdd- 

above together with the typical Lorentz factors result in the jet 



• The values of Mj Ct 



powers Pj c t ~ 10 — 10 erg s . 

• Given the bolometric luminosities Lboi, we have (Pj e t/^boi) ~ 10. 
Taking the values above, we can roughly estimate the fraction of accreted mass that is 

channeled into the jets in LINERS, (Mj et /M) ~ 10~ 3 . Similarly, we can estimate the 



efficiency of jet production as (P- ]ct / (Mc 2 )) ~ 1%, in rough agreement with other esti- 



mates for LLAGNs (jAllen et al.l 2006). These values provide useful indicators of the rele- 
vant feeding and feedback properties of LINERs, and by extension of the whole LLAGN 
population; this is particularly relevant in the light of this symposium. 

Finally, our work allows us to draw a link between the supermassive black holes in LIN- 
ERs and the quiescent b east in our galaxy, Sgr A*. By mod elling the SED of Sgr A* with 
current ADAF models, Yuan. Quataert fc Naravanl ( 20031 ) estimated M ~ 10~ 5 MEdd- 
This value is two orders of magnitude below the typical accretion rate of LINERs that 
we estimated. Therefore, we could say that if Sgr A* accreted at a rate 100 times its 
present rate, it would "light up" and presumably become a LINER. 



5. Conclusions 

We were able to successfully model the SEDs of 24 LINERs in the context of a coupled 
ADAF-jet scenario. While the the radio emission is dominated by the relativistic jet, the 
ADAF dominates in the band 1 mm - 100 /xm and the X-ray radiation can be dominated 
by either the ADAF, the jet or a combination of both. 
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We find that strong jets are implied by our modelling, for which the kinetic power 
considerably exceeds the radiated power ((Pj Ct /Lboi) ~ 10). Furthermore, we obtained 
estimates of the fundamental parameters of the central engines of LLAGNs which can be 
useful in studies of the feeding/feedback properties of AGNs, such as the mass accretion 
rates, jet powers and mass-loss rates in the jet. 

Finally, we would like to point out that our SED models provide a library of tem- 
plates for the compact emission of LINERs, which can be useful for studies of e.g. the 
nuclear emission of stellar populations, dust and PAH features. A detailed description 
of the results above and the SED models will be available in Nemmen et al. (2010), in 
preparation. 
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